High-resolution Spectroscopy of [Ne II] Emission from TW Hya 
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ABSTRACT 

We present high-resolution echelle spectra of [Ne II] 12.81/im emission from 
the classical T Tauri star (CTTS) TW Hya obtained with MICHELLE on Gem- 
ini North. The line is centered at the stellar radial velocity and has an intrinsic 
FWHM of 21 ± 4 km s~^. The line width is broader than other narrow emission 
lines typically associated with the disk around TW Hya. If formed in a disk, 
the line broadening could result from turbulence in a warm disk atmosphere, 
Keplerian rotation at an average distance of 0.1 AU from the star, or a photoe- 
vaporative flow from the optically-thin region of the disk. We place upper limits 
on the [Ne II] emission flux from the CTTSs DP Tau and BP Tau. 

Subject headings: stars: pre-main sequence — stars: planetary systems: proto- 
planetary disks matter — infrared: stars: TW Hya, DP Tau, BP Tau 
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1. INTRODUCTION 

X-ray and UV emission from classical T Tauri stars (CTTSs) can affect the chemical 
structure of their protoplanetary disks by ionizir ig the disk surface , dissociating molecules 



and heating the g as and grains (see reviews by iNajita et al.l 120071 : iDuUemond et al.l 12007 



Bergin et al.ll2007l ). Direct evidence of irradiated disks is seen in far-UV H2 emission, which 



requires excitation by H I Lyg irradiation and can be attributed to emission from the inner 



disk of some CTTSs (IHerczeg et al.ll2006l ). Other observational evidence for FUV- and X- 
ray-irradiated disks around CTTSs has been more circumstantial. The underabundance 



Bergin et al. 


2003 


)• 


Barv et al. ( 


2003) 



disk surface enough to produce emission in the H2 1-0 S(l) line. The excess FUV continuum 
emission may result from the electronic cascade and dissociation following excitation of H2 
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(Herczeff et al. 


2004; 


Berein et al. 


2004) 



can produce a photoevaporative wind, which could be visible in the l ow- velocity component 



of some optical forbidden lines (see review by lDuUemond et al.l 120071 ) 



Recently, iPascucci et al.l (120071 ) and iLahuis et al.l (120071 ) may have found a new probe 
of gas at the disk surface by detecting emission in the [Ne II] 12.81/im fine-structure line in 
low -resolution Spitzer IR S spectra of CTTSs. [Ne II] was dete cted in four o f six t argets in 
the iPascucci et al.l (120071 ) sample and in 15 of 76 targets in the llahuis et al.l J2007h sample. 
Relative to the ~ 1 Myr old CTTSs in Taurus, the six stars observed by IPascucci et al. 



( 120071 ) have faint mid-IR continuum emission, are X-ray bright, and have mass accretion 
rates an ord e r of m agnitude lower than the average CTTS of similar spectral type. The 
Lahuis et al.l (|2007l ) sample and detections are comprised mostly of T Tauri stars with an 
IR excess. The only class I source in their sample was also detected in [Ne II], and several 
Herbig Ae/Be s tars were undetecte d in [Ne II]. The sources both with and without [ Ne II ] 
detections from lLahuis et al.l (120071 ) show a range of X-ray luminosities. iLahuis et al.l (120071 ) 
also detected [Ne III] emission from one source (Szl02), and four of the six targets with H2 
0-0 S(2) emission also show [Ne II] emission. 



Pascucci et al.l ( 120071 ) and ILahuis et al.l (120071 ) both argue that the irradiated disk sur- 
face is the most likely origin for thi s emission. Tw o ioniz ation paths can ionize Ne and 



thereby produce the [Ne II] emission. iGlassgold et al.l (120071 ) propose that the line is formed 
because of Ne ionization by K-shell absorption of stellar X-rays at energies > 0.9 keV. The 
observed [Ne II] fluxes and the [Ne III]/[Ne II] flux ratio and lower limits are consistent 
with predictions from this X-ray ionization model. Alternatively, EUV photons shortward 
of the Ne I ionization edge at 575 A can produce [Ne II] emission if the EUV emission is 
able to penetrate through neutral gas in the accretion column and wind. If the disk origin 
is confirmed, then the [Ne II] 12.81/im line would be a direct probe of the irradiation of the 
protoplanetary disk by the central star. 

In this paper, we present high-resolution echelle spectra of the 12.81yum region of the 
CTTSs TW Hya, BP Tau, and DP Tau obtained with MICHELLE on Gemini North. TW 
Hya is a nearby 10 Myr old CTTS with a disk that is believed to be optical ly-thin within 4 AU 



from the star based on models of its infrared spectral energy dis tribution JCa 



vet et al.ll2002l) 



2007h. BP Tau 



and an inner hole and bright ring imaged at 7mm with the VLA (iHughes et al 
and DP Tau are 1 Myr old CTTSs with optically-thick inner disks (IKenyon fc Hartmann 
I995I : iMuzeroUe et al.ll2003l ). All three stars are X-ray bright (Table 1). We find strong 
[Ne II] emission from TW Hya and place upper limits on the [Ne II] flux from DP Tau and 
BP Tau. We discuss the origin of the emission in the context of profiles of other emission 
lines from TW Hya and suggest that the line is formed in a disk, although formation in a 
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heated photosphere cannot be ruled out. If the [Ne II] emission is formed in a disk, the line 
width could result from Keplerian rotation at a distance of 0.1 AU from the star, tran-sonic 
turbulence, or detecting the front and back of a photoevaporative flow. 



2. OBSERVATIONS AND DATA REDUCTION 

We used the echelle mode of MICHELLE on Gemini North to observe the CTTSs TW 
Hya, DP Tau, and BP Ta u on 9 March 200 7. MICHELLE has a 320x240 Si:As array 



with a plate scale of 0'.'201 (IGlasse et al.l 119971 ). When used with the echelle spectrograph 



MICHELLE has a spectral resolution of ~ 30, 000. We centered the targets on the slit using 
images obtained with the Si-5 filter, which has an effective wavelength of 11.7 yum. We used 
a 2-pixel wide slit and the Si-6 filter to obtain spectra centered at 12.815/im and covering 
~ 0.049 iim. We nodded 10" off the slit at 35s intervals. Our total on-source integration was 
40.4 min for TW Hya, 20 min for BP Tau and 22.5 min for DP Tau. These three targets 
were observed at median airmasses of 1.8, 1.3, and 1.9, respectively. Observations of a CMa 
and the asteroid Iris were used as telluric standards. The seeing varied from 0'.'6 to 0'.'85, as 
measured in the cross-dispersion spatial profile from our unresolved telluric standards. 

In each set of four nod positions, the counts in each column of the off-source images 
are scaled to and then subtracted from the counts in the same column of the on-source sky 
background. This scaling is necessary because the second derivative of the sky brightness 
versus time was non-zero. We correct for fringing in the cross-dispersion direction by fitting a 
high-order polynomial to the residual background in each column. We extract the spectrum 
from a 7-pixel window centered on the source, and then subtract background sky spectra 
obtained from 20-pixel windows on both sides of the source. Fringing in the spectral direction 
is not detected or corrected for because of the low S/N in the continuum (Table 1). 

We estimate the telluric correction for each observation by interpolating the standard 
spectra of Iris and a CMa at intermediate airmasses. The atmospheric correction is applied 
to each nod prior to combining the images. Several O3, CO2, and H2O lines are present in 
the wavelength range of our spectra. We calibrate wavelengths to an accuracy o f ~ 2 km s~ ^ 



per pixel using a 2nd-order fit to sky emission lines identified with ATRAN (ILordl Il992l ). 
Each pixel covers ~ 3.6 km s~^. We estimate a resolution of 3.2 ± 0.5 pixels (11.5 ± 1.8 
km s~^) from the FWHM of four weak, narrow sky emission lines. 

TW Hya was also observed on 12 April 2 006 during re-commissioning of the East Arm 



Echelle Spectrograph (ILibbrecht Sz Peril 119951 ) on the 5m Hale Telescope at Palomar. The 



instrument is fed by a bundle of fibers that together subtend 1'.'5 and illuminate a 2048 x 2048 
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pixel CCD. The spectrum covers most of the 4000-9000 A range at a resolution of 15,000- 
45,000, depending on location on the chip. The relative wavelength solution in each order 
was obtained using a ThAr lamp. The absolute wavelength solution is determined separately 
for each echelle order by matching photospheric lines to the radial velocity of TW Hya and 
is accurate to ~ 3 km s~^. 

No sufficient flat field was obtained during the commissioning. However, observations 
of the white dwarf Feige 34 and the AO star HD74198 obtained in the same night yield 
spectra with a count rate that peaks near the center of each order and falls off smoothly 
to the edge of the detector. The absence of flat-fielding has a negligible effect on the line 
profile shapes and equivalent width measurements because they are extracted from narrow 
spectral regions. We use observations of several other stars to identify the approximate 
flux in telluric emission lines, which are located at -20 km s~^ relative to the photospheric 
velocity of TW Hya. The telluric and stellar components of the [O I] lines at 5577, 6300, 
and 6363 A can therefore be separated. The Na I D emission lines are undetected, with 
a sensitivity limited by strong telluric emission. We calculate line fluxes from equivalent 
widths by estir nating the continuum flux level from non - simultaneous archival HST/STIS 
optical spectra (IHerczeg et al.ll2004j : iRoberge et al.ll2005l ). 



The heliocentric radial velocity of TW Hya is 12.2 ± 0.5 km s~^ (Weintraub et al. 
2000; see also Kastri er et al. 1997, de la Reza et al. 1989) and of BP Tau is 15.8 km s~^ 
( iHerbig fc Belli Il995l ). The heliocentric radial velocity of DP Tau has not been measured. 
We estimate a radial velocity of 17 ± 5 km s~^ by aver aging the radial velocities of the T 
Tauri stars GO Tau and V955 Tau (IHerbig fc Belllll995l ). which are both about 10' from DP 
Tau. 



3. CHARACTERIZING THE [Ne II] EMISSION 



Figure 1 shows our MICHELLE spectra of TW Hya, BP Tau, and DP Tau. The 
continuum is detected from each target. We measure e quivalent widths for the [Ne II] flne- 



structure line (^Pg/s - at 12.81355 ±0.00002 fim (lYamada et al.lll985h of 62 ± 11 A for 



TW Hya and 2a upper limits of < 10 and < 11 A for non-detections from BP Tau and DP 
Tau, respectively. The upper limits are calculated by assuming the same line proflle as that 
of TW Hya, and would be larger if the hypothetical lines from these objects were broader. 
We convert equivalent widths to (5.9 ± 1.1) x 10"^^, < 0.7 x 10"^^, and < 1. 1 x 10"^"^ erg 



cm 



^ using 12.8jUm c ontinuum flux measurements from Spitzer IRS spectra (jUchida et al. 



2OO4J : iFurlan et al.ll2006l ). These fluxes assume that the mid-IR continuum and the [Ne II] 
emission are not spatially extended beyond the seeing. 
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The [Ne II] line profile from TW Hya is consistent with a single Gaussian profile centered 
at — 2±3 km s~^ in the stellocentric frame with a FWHM of 6.7±0.8 pixels (24±3 km s~^). 
Deconvolving the line width from the instrumental resolution yields an intrinsic line FWHM 
of 21 ± 4 km s~^. Two components (a broad component at the stellocentric velocity with 
a flux enhancement on the blue side of the line) are visually apparent in the line profile. 
However, based on the S/N detected in the continuum (Table 1) the profile shape is not a 
significant deviation from a single Gaussian profile. The telluric CO2 line did not divide out 
sufficiently for some of the TW Hya observations and results in increased uncertainty on the 
blue side of the line profile. The [Ne II] emission is more compact than the 0'.'75 angular 
resolution (40 AU at the distance of TW Hya), which was estimated from both the mid-IR 
continuum of TW Hya and the telluric standard a CMa, which was observed shortly before 
TW Hya. 

Table 2 and Figure 2 compare the [Ne II] emission to a sample of other emission lines 
from TW Hya, excluding those significantly affected by wind absorption (e.g., H I A6563, O I 
A1302, C II A1335, He I A10830, see Herczeg et al. 2002, Edwards et al. 2006). These lines 
show a range of emission profile shapes, and were observed at resolutions similar to that of our 
MICHELLE spectra. The velocity and FWHM listed in Table 2 are measured from Gaussian 
fits to the emission lines. The energy of the upper level (i?') for these transitions (0.60-23.1 
eV) are all substantially larger than that of the fine-structure [Ne II] line (0.10 eV). Most 
lines, including the Ca II IR triplet and the [O I] A6300, 6363, show narrow emission centered 
at the stellocentric velocity. No high-velocity (> 100 k m s~^) or low-velocity (< 20 km s~^ ) 



component of other optical forbidden lines are detected (lHamannlll994l : lHartigan et al.lll995l ). 
Several lines show broad redshifted emission profiles, presumably produced by the accretion 
flow. Many of the near-UV Fe II lines have both narrow and broad components centered at 
the stellocentric velocity. 



4. DISCUSSION 

4.1. Comparison To Other Emission Lines Prom TW Hya 

Emission lines from CTTSs may be formed at or near the disk surface, in the accretion 
funnel flow and shock, in outflows, in the stellar chromosphere, or in the accretion-heated part 
of the stellar photosphere. We first rule out formation of the [Ne II] emission in the accretion 



A1549, N V A1240, and VI A1035 doublets ( 


A 


encar & Batalha 


2002; 


Herczee et al. 


2002; 


Lamzin et al. 


2004: Johns-KruU & Herczee 


2007 


Yang et al.l 


2007 


), have much broader line 



widths than the [Ne II] emission from TW Hya and are redshifted from the stellocentric 
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velocity. Second, we rule out fast outflows because the [Ne II] line is centered at or near 
the stellocentric velocity and is limited to low velocities. Formation in a very slow outflow 
may be possible if both the front and back sides of the flow are visible. Third, we rule out a 
chromospheric origin. Some lines from Table 2 are attributed t o chrornosph eric emission and 
seen from main-sequence dwarfs and weak-line T Tauri stars (jAyresll2005l ). but these lines 
are narrower than the [Ne II] line. Several other emission lines, including the C I A1656 and 
C II] A2326 multiplets, CI I A1351.7, and some NUV Fe II lines, have slightly larger widths 
than the [Ne II] emission and may form in a heated stellar photosphere. However, these lines 
have larger energy upper levels than the flne-structure [Ne II] line and are unlikely to trace 
the same gas. More likely is the possibility that the [Ne II] line forms either in the inner 
disk or in a photoevaporative flow from the optically thin region of the disk, as suggested by 
the excitation energy (0.10 eV), width (21 km s~^), and velocity centroid (—2 ± 3 km s~^) 
of the line. 

Collisional ionization equilibrium would require temperatures of > 30, 000 K to produce 
[Ne II] emission. All identifled lines from TW Hya that probe > 20, 000 K have a larger 
width (> 40 km s~^) than the [Ne II] line. We therefore infer that the [Ne II] emission 
is not likely coUisionally ionized. Instead, the emission is likely ionized by X-ray or EUV- 
irradiation of gas at or near the disk surface. Indeed, Glassgold et al. (2007) speculated 
that X-ray irradiated disks would produce strong [Ne II] emission. The observed width of 
the [Ne II] line is larger than that shown by Glassgold et al. (in their Fig. 5). However, 
Glassgold et al. (2007) considered only the emission that would arise from a non-turbulent 
X-ray irradiated disk that extends beyond 1 AU. 

Irradiated disks are believed to produce the molecular line emis sion detected from TW 



Hya in the FUV H ^ , near-IR Ho, and fundamental CO t ransitions (IWeintraub et al.l 12000 



Herczeg et all l2002l : iRettig et al.l 12004 : ISalvk et all l2007l : iNaiita et al.l l2007h . The emission 
arises from warm (500 — 3000 K) gas, and the FUV H2 lines in particular require irradiation 
by strong FUV emission from the central star. The line widths are > 16, 8, and 14 km s~^ 
narrower for H2 1-0 S(l), FUV H2, and fundamental CO lines, respectively, than the [Ne II] 
line. The narrow width of the H2 1-0 S(l) suggests formation at large radii. 

Interestingly, the widths of the fundamental CO, FUV H2, and [Ne II] lines increase with 
their approximate formation temperatures of ~ 500 — 1500 K (Salyk et al. 2007; Najita et al. 
2003), ~ 2500 K (as observed for TW Hya; Herczeg et al. 2004), and 4000-10,000 K (based 
on theoretical modeling; Glassgold et al. 2007; HoUenbach et al. 2000), respectively. This 
trend might indicate a common origin for the lines in a disk with an outwardly decreasing 
temperature gradient (e.g., Carr et al. 2004). If the [Ne II] line width results from Keplerian 
disk rotation, then the emission arises at an average distance of ~ 0.1 AU of the star. 
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assuming a 7 ± 1° inclination of the disk around TW Hya (iQi et al.ll2004j ). The trend might 
also indicate a common origin for the lines in a turbulent disk atmosphere with a vertically 
increasing temperature gradient,- i.e., increasing (transonic) turbulent broadening higher in 
the disk atmosphere. For example, the line width of the [Ne II] emission could be explained 
entirely by turbulent broadening in 10,000 K gas. With a combination of Keplerian rotation 
and turbulence, the line could form at disk radii larger than ~ 0.1 AU and in gas cooler than 
10, 000 K. Alternatively, the width of the [Ne II] line could also be explained if the emission 
arises from a photoevaporative flow from the optically thin region of the disk (at < 4 AU). 
We explore these possibilities using simple scaling arguments to point out the challenges 
associated with these interpretations. 



4.2. X-ray Excitation of the Inner Disk 



The [Ne II] flux is converted to a mass from the Einstein A-value of 8.59 x 10~'^ s 
(iKramida fc Navdl2006l ). a solar Ne/H abundance ratio of 8.9 x 10~^ ( lAnders fc Gervesse 
19891 ). a Ne ionization fraction of C, and an average atomic mass per atom of 1.4. In LTE the 
fraction of Ne II in the ^^3/2 level is 0.14 and 0.28 at 1000 and 5000 K if the electron density 
is larger than the critical density, 10^ cm~'^, as assumed here. If the line is optically-thin and 
not attenuated by dust within the disk, then the gas mass is 0.5 — 1.5 x 10~^C~^ where 
this range includes uncertainties in the line flux and the population of the upper level. If we 
assume that the [Ne II] emission arises in the annulus of an empty cylinder with a height H, 
radius i?, and width dR < R, then the total hydrogen column density of the layer measured 
parallel to the disk is N{}i) ~ 1.5 x 1022^^0^^- 



-1 —2 

cm . 



Either the disk atmosphere or the photoevaporative flow scenario requires that the ion- 
izing photons penetrate through the emission region, which constrains the neutral hydrogen 
column density, A^(H I), as r aeasured radially fro m the central star. Based on the Chandra 
X-ray spectrum of TW Hya (IKastner et al.ll2002l ). 23% of photons with energies > 0.9 keV 
penetrate to depths of A^(H I)=10 
penetrate to A^(H I)=10 
A are attenuated for A^(H I)> 10^° and 10^^ 



cm~^, but only the hardest 3% of those photons can 
cm~^. However, almost all EUV photons longward of 100 and 500 
cm~^, respectively. 

If we set R < 0.2 AU, then A^(H)~ lO^^C"^ cm~^. Even if the ionization fraction is 
high, A^(H I) will be large enough to require the ionization of Ne by X-rays and a disk height 
comparable to the disk radius. Such a geometrically thick disk may also be required to e xplain 
the strength of the H2 fluorescence in FUV spectra of TW Hya (IHerczeg et al.l 120041 ). The 
scale height of gas in h ydrostatic equilibriu m at i? = 0.2 AU from TW Hya can be estimated 
as 0.04 (yp^) AU. Eisner et al. ( 2006 ) calculated that the dust temperature is ~ 1100 
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K at the inner dust truncation radius (~ 0.06 AU). However, si nce only a sraall am ount of 
micron-sized dust is present within the 4 AU dust clearing radius ( ICalvet et al.ll2002l ). the gas 
temperature at these radii may be significantly hotter than the dust temperature. Models of 
optically-thick disks with viscosity parameter a = 0.01 su ggest a gas surface density of 25-90 
g cm-2 at 0.1 AU from TW Hya ( 



D'Alessio et al.lll998l ). The mass of the [Ne II] emitting 



region would be about ^% of the total disk mass within 0.2 AU of the star. The [Ne II] 
emission in this scenario traces a thin surface layer of the disk higher than one scale height 
above the disk midplane. This interpretation of [Ne II] emission is challenging if ( < 0.1 
because the larger total mass traced by the [Ne II] emission would require an even larger 
geometrical disk height. 

This situation can be ameliorated if turbulent broadening contributes to or dominates 
the [Ne II] line width. If the disk atmosphere is turbulent and the turbulent velocity is close 
to the sound speed, then the FWHM of a line would be ~ 21 km s~^ at 10, 000 K, consistent 
with the observed line width. At this temperature the [Ne II] emission could arise from a 
wide range of radii, perhaps extending out to 10-20 AU. iGlassgold et al.l (120071 ) suggest that 
at 1-10 AU, the gas within N{H) = 10^° cm"^ of the disk surface will be heated to 4000 K. 
At this temperature the average radii for formation of [Ne II] emission would be ~ 0.5 AU 
from the star. 

Approximately transonic turbulence in disk atmospheres is inferred from the intrinsic 
line broadening of CO overtone emission from accretion disks (Carr et al. 2004; Najita et al. 
1996). The dissipation of turbulence in disk atmospheres could also produce the extra heat- 
i ng that is needed to explain the strength of CO fundamental emission from T Tauri disks 
( iGlassgold et al.ll2004l ). While these observations refer to emission from cool, molecular disk 
gas that is located deeper in the disk atmosphere, similar transonic turbulence may charac- 
terize warmer, atomic gas higher up in the disk atmosphere where the [Ne II] emission may 
arise. Indeed, simulations of weakly magnetized disks find that MHD turbulence generates a 
warm 10, 000 K) disk coro na with a transonic velocity dispersion above several disk scale 
heights jMiller fc Stonelboooh . 



4.3. Photoevaporation of an EUV-irradiated Disk Surface 



On the other hand, in the disk dissipation models of Alexander et al. (2006ab) EUV irra- 
dation of the disk produces a 10, 000 K fiow with a velocity of ~ 10 km s~^ from the disk sur- 
face, near the gravitational radius Vg ( see also HoUen bach & Gorti 2007). Analytical formula- 

AURi 1 AU 



tions for assuming an abiabatic gas (jLiffmanll2003l ) suggest that Vg f» 1.4 
for a CTTS with the mass of TW Hya. In the photoevaporation simulations by iFont et al. 



Mq' 
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( 120041 ). most of the mass loss occurs inside ~ 6 AU for a source l ike TW Hya. While these 
expectations are for a radially continuous gas disk, ICalvet et al.l (120021 ) find that at 10 fim 
the disk is optically thin within 4 AU of TW Hya. Any additional frontal illumination of 
the inner edge of the optically-thick disk at 4 AU by FUV and X-ray photons may provide 
an additional heating source that could help puff up the disk or help launch the photoevap- 
orative flow. If the [Ne II] emission is formed in a photoevaporative flow within 4 AU of the 
star, the apparent line broadening could result from detecting [Ne II] emission from both 
the front and back sides of the disk. The proflle could be slightly blueshifted if some [Ne II] 
emission occurs in a photoevaporative flow beyond th e dust clearing radius, w here emission 
on the back side of the disk would not be detectable. [Alexander et al.l (l2006al ) predict that 
the photoevaporative flow will have velocities > 10 km s~^ at distances greater than 4 AU 
from the star, implying the possibility of a stronger blueshifted asymmetry than is observed. 



^ 4 AU, then the height of the 

1.5xl020cm-2 If C = 1 and 



If we assume the photoevaporative flow is produced at R 
flow (including both front and back sides of the disk) is H — (n{r) 
A^(H)~ 10^^ cm~^, then H = 15 AU. The v elocity of a photoevapor ative flow is expected to 
be f ~ 10 km s~^ at 4 AU above the disk (lAlexander et al.ll2006al ). The crossing time for 



the evaporating gas across the emission region on either side of the disk {H/2) is then 3.5 
yr. The estimated mass loss rate is ~ 3 x 10~^°/ Mq yr~^, where / is the fraction of 
[Ne II] emission that arises in the photoevaporative flow. 

Such a low-velocity (~ 10 km s^^) photoevaporative flow would be distinct from the 
observed high-velocity (—50 to —200 km s~^) wind that dominates blueshifted absorption 
proflles of atomic and singly-ionized lines from TW Hya (Edwards et al. 2006; Johns-KruU 
et al. 2007). The estimated photoevaporation mass loss rate of about one-flfth the stellar 
mass accretion rate is similar to the ~ 3 x lO^^'' Mq yr~^ estimate for a radially-continuou s 
disk that is irradiated with an EUV photon flux of ~ 10^^ phot s~^ (iHoUenbach et al.lll994j ). 
Several parameters in this simple scaling argument, including /, are uncertain and together 
could suppress our estimate for the mass loss rate. This scenario constrains the iV(H) in 
the photoevaporative flow to a narrow range. A smaller A^(H) would increase the volume 
and decrease the Ue below the critical density by increasing H, where any flow is likely less 
collimated. A larger A^(H) decreases the height and crossing time, thereby increasing the 
estimate of the mass loss rate. 

The photoevaporation scenario is contingent upon some EUV emission reaching the disk 
surface. The ph oton flux at < 912 A fr om TW Hya of ~ 10^^ — 10^^ phot s~^ is dominated 
by line emission (lAlexander et al.ll2005l ). Since both the ~ 10^ K FUV l ines and 10^ K X-ra y 



lines are attributed to the accretion flow (jJohns-KruU &: Herczegl 120071 : iKastner et al.ll2002l ). 
this estimated EUV emission likely applies to emission produced near the accretion shock. 
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Whether this EUV emission reaches the disk depends on a small neutral hydrogen column 
density. A smaller amount of EUV emission is expected to be produced by the transition 
region in magnetic structures and may reach the disk, de pending on the geometry of the 
accretion flow. Based on an analysis of FUV H2 emission, iHerczeg et al.l (120041 ) suggested 
that an absorbing column of A^(H I)~ 5 x 10^^ cm~^ H I, perhaps in a stellar outflow, is 
present between the H I Lya emission from TW Hya and the warm H2 in the disk. If 
A^(H I)< 2 X 10^*^ between the EUV and [Ne II] emission resions, then 5-15% of the EUV 
photons could reach the disk and lead to some photoevaporation. A larger intervening A^(H I) 
would prevent EUV photoevaporation. 



4.4. Conclusions 



Our lower luminosity limits for [Ne II] emission from BP Tau and DP Tau are com- 
parable to the line luminosity from TW Hya, despite having mass accretion rates an order 
of magnitude larger than that onto TW H ya. [Ne II] is de t ected from several CTTSs with 
low mass accretiori rates both here and in iPascucci et al.l (120071 ). Of the few stars in the 
Lahuis et al.l (120071 ) sample with measured mass accretion rates, [Ne II] emission is detected 
from several weak accretors but not from several strong accretors. If weak accretors do show 
stronger [Ne II] emission than strong accretors, this effect may be explained because the 
intervening iV(H I) in the accretion flow and wind should scale with the mass accretion and 
mass loss rates. However, a similar effect may be seen for X-ray-irradiation models if older 
CTTSs are more X-ray luminous than younger CTTSs. 

The [O I] emission from TW Hya may be produced by recombination subsequent to 
photoionization. T he [O I] lines are centered at the stellocentric velocity, do not correlate 
with optical veiling (lAlencar fc Batalhall2002l ). and are not seen from WTTSs. [O I] enaission 



is thought to ari se from the disk surface of some Herbig Ae/Be stars (IBohm fc Catalalll994 



A eke et al.ll2005l ). [O I] emission produced subsequent to OH dissocia tion by FUV emission 
from the central star will be undetectable from disks around CTTSs (IStorzer fc HoUenbach 
2OOOI ). However, models of EUV or X-ra y irradiated disks may explain some [O I] emission in 
the low- velocity components of CTTSs (IFont et al.ll2004j ). That the [O I] lines are narrower 
than the [Ne II] hne suggests that, if the [O I] emission is produced in a disk, it traces less 
turbulent gas or gas at larger disk radius than the [Ne II] emission. 

High-resolution spectroscopy of [Ne II] emission from TW Hya suggests formation in a 
disk irradiated by X-ray or EUV photons, although formation in a heated photosphere cannot 
be ruled out. The origin of the line broadening and the ionization path can be identified 
by high-resolution spectra of additional sources. If the line broadening is dominated by 
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Keplerian rotation, the [Ne II] line will appear much broader from stars with disks viewed 
at higher inclinations. If the line broadening is dominated by turbulence, the observed line 
width from disks viewed at higher inclinations will only be modestly larger. If formed in a 
photocvaporative wind, the [Ne II] emission from stars with optically-thick inner disks will 
likely appear blueshifted, as only the front side of the disk will be detectable. The description 
of [Ne II] emission from TW Hya should apply generally to stars with small mass accretion 
rates. [Ne II] emission from stars with large mass accretion rates may also be produced in 
shocks associated with outflows. 
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Table 1: Source Properties 



Star 


SpT 


d 


Vr 


Lx 


M* 


R, 


M 


-'^12.8 


S/N'' 


EW ([Neil]) 


LncII 


Refs 






pc 


km s^^ 


1O-4L0 


Mq 


Rq 


IQ-^Mq yr-i 






A 


io-bl© 




TW Hya 


K7 


51 


12.2 


1.4^= 


0.7 


1.0 


1.8,0.5°' 


0.95 


1.3 


62± 11 


4.8 


1,2,3,4 


BP Tau 


K7 


140 


15.8 


1.8"= 


0.5 


2.0 


24 


0.7 


1.5 


< 10 


< 4 


5,6,7,8 


DP Tau 


MO 


140 


17f 


0.8<= 


0.5 


1.4 


9.0 


1.0 


1.8 


< 11 


< 7 


5,6,7,8 



^) at 12.8/im from Spitzer IRS spectra 
''S/N per pixel in the continuum in our MICHELLE observations 

''0.45-2.25 eV luminosity from high-resolution XMM spectra (Stelzer fc Schmitt 2004 ). 

'^1.8 from Alencar fc Batalha (2002) and Herczee et al. (2004), adjusted for 51 pc distance. 0.5 from MuzeroUe et al. ('2000) 
'=0.3— 10 eV from models based on low-resolution XMM spectra (Briggs_et_aL ^007) . 
■''Average radial velocities of GO Tau and ■V955 Tau from Herbig fc Bell (1995). 

IMamaiek (2005) 2Uchida et al. (2004) 3Weintraub et al. (2000) 4'Webb et al. (1999) 

SValenti et al. (1993) 6Herbig &: Bell (1995) 7Briggs et al. (2007) SHartmann et al. (1998) 
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Table 2: Selected Emission Lines from TW Hya, in order of FWHM"-'' 



Line 


^oba 


E' 


R 


V 


FWHM'' 


EW 


Flux'^ 


Formation'* 


ID 


A 


cV 


A/AA 


km s^-"- 


km s~^ 


A 






Hi 


21218 


0.60 


49000 





< 6 


0.02 


1 


D 


CO (fundamental)-^ 




24000 





8.3 


- 


- 


D 


Si I 


2882.551 


5.1 


30000 


1 


11 


0.14 


0.96 


P 


[O I]" 


(5577.571) 


2.0 


36000 


(0) 


(12) 


(0.11) 


< 1.4 


- 


[Ol] 


6300.549 


2.0 


34000 


-1 


12 


0.67 


8.7 


D? 


[O I] 


6364.050 


2.0 


26000 





13 


0.23 


3.5 


D? 


[S 11]" 
[N 11]" 


(6730.81) 


1.8 


20500. 


(0) 


(12) 


< 0.06 


< 0.9 


- 


(6583.46) 


1.9 


21000 


(0) 


(12) 


< 0.07 


< 1.0 


- 


[Fe II]" 


(8616.952) 


1.7 


27500 


(0) 


(12) 


< 0.07 


< 1.2 


- 


Na Is 


(5889.951) 


2.1 


22000 


(0) 


(12) 


< 0.45 


< 4.5 


- 


Na If 


(5895.924) 


2.1 


22000 


(0) 


(12) 


< 0.45 


< 4.5 


- 




FUV'' 


1-2' 


25000 





14 


- 


- 


D 


Fell 


2371.319 


5.6 


30000 


-1 


15 


0.17 


0.76 


? 


Fell 


2665.569 


8.0 


30000 





15 


0.09 


0.55 


7 


Ca II 


8542.362 


3.2 


16000 


-2 


15 


1.2 


21 


P+A? 


Ca II 


8662.464 


3.1 


43000 


-1 


16 


0.84 


17 


P+A? 


Mg II 


2791.722 


8.9 


30000 


1 


16 


0.27 


0.19 


P? 


[Ne II] 


128140 


0.10 


31000 


-2 


22.5 


62 


5.9 


D 


Fe II 


2626.565 


4.8 


30000 


1 


25 


1.3 


6.3 


P/H? 


{2y 


2626.557 


4.8 


30000 





130 


2.1 


10.3 


H? 


CI I 


1351.715 


9.3 


30000 





41 


0.18 


0.27 


H? 


Cl 


1656.346 


7.5 


30000 


1 


41 


1.5 


4.6 


H? 


ClI] 
He P 


2326.231 


5.3 


30000 


1 


42 


3.8 


19 


H? 


5876.060 


23.1 


19000 


10 


50 


0.81 


12 


A 


Si 


1296.231 


9.6 


30000 





52 


0.95 


1.3 


A? 


nI] 


1411.966 


12.4 


30000 


-7 


83 


1.9 


2.8 


A? 


He IP 


1640.601 


48 


30000 


20 


88 


310 


94 


A 


OlII] 


1666.093 


7.5 


30000 


-23 


113 


0.8 


1.6 


W7 


cm 


2297.842 


18.1 


30000 


22 


127 


0.8 


3.2 


A 


H7 


4341.214 


13.1 


32000 


39 


195 


8.1 


81 


A 


H/3"' 


4861.713 


12.8 


32000 


11 


199 


34 


44 


A 


H<5 


4102.464 


13.2 


32000 


41 


201 


2.9 


23 


A 


cm 


1548.757 


8.0 


30000 


110 


300 


93 


186 


A 


N V 


1239.068 


10.0 


30000 


60 


340 


15 


30 


A 



"Intrinsic FWHM after deconvolution from the instrumental resolution. 
"Assumed properties for undetected lines in parentheses. 

c-|^g— 14 gj,g (.ijj-2 g-l 

■'W=Wind, P=photosphere, H=Heated Photosphere, A=Accretion funnel, D=Disk 

=From Weintraub et al. (2000). 

■'■ Averaged CO FWHM from Salyk et al. (2007). 

f Not detected here. ^lencar_&_^atalha (2002) found narrow Na I D lines at the stelloccntric velocity. 
''Coadded FUV H2 profiles by coadding many of the strongest lines from Herczeg et al. (2002). 
'Excited from many levels between 1—2 eV (Herczeg et al. 2002). 
^NUV Fe II lines from TW Hya have narrow and broad components. 

^Properties are for bright, narrower component. Excess weak redshifted emission is also present. 

'Line center and width calculate from a Gaussian fit to a non-Gaussian profile, see Johns-KruU & Herczeg (2007). 

'"Corrected for optically-thin wind absorption 
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Fig. 1.— The MICHELLE spectra of TW Hya, BP Tau, and DP Tau, shifted to the ap- 
propriate stellar velocity. The [Ne II] line is strong from TW Hya but not detected from 
BP Tau or DP Tau. The mid-IR continuum is detected from all three sources. The bottom 
panel shows the spectrum of the telluric standard a CMa shifted to the observed velocity of 
TW Hya. 
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Fig. 2. — A comparison of the [Ne II] line profile with many optical and UV lines from TW 
Hya. The [Ne II] line is much more narrow than lines associated with the accretion funnel 
(including the He I A5876 line shown here), but is broader than lines formed in the disk (H2 
and possibly [O I] emission shown here) and in the photosphere (Si I emission shown here). 
Telluric [O I] emission is identified by the Earth symbol. 



